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ABSTRACT 

The  emplacement  and  recovery  of  large  deep  ocean  cable  systems 
containing  in-line  packages  require  a  knowledge  of  the  dynamic  tension 
and  motions  of  the  system  under  the  influence  of  surface  ship  motions, 
subsurface  ocean  currents  and  other  external  loads.  A  two-dimensional 
lumped  mass  model  was  developed  to  simulate  simply  connected  cables  and 
in-line  packages.  Cable  tensions  and  mass  poin  position  and  velocity 
are  determined  by  a  finite  difference  method  using  a  predictor  -  corrector 
technique.  The  resulting  computer  program  is  applicable  to  single  pay- 
load  lowering,  cable  laying,  deep  sea  mooring  and  deployment  of  large 
undersea  cable  structures. 
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INTRODUCTION 


Ocean  mooring  and  lifting  are  necessary  surface  support  capabilities 
for  the  emplacement  of  seafloor  structures  and  underwater  cable  systems. 
To  prope, ’v  select  mooring  and  lifting  lines,  predictions  of  line  ten¬ 
sions  throughout  the  various  operational  phases  are  required.  The 
selection  of  undersi ged  lines  due  to  lack  of  dynamic  stress  information 
can  result  in  catastrophic  failures.  Therefore,  methods  of  analyzing 
cable  systems  are  needed  for  proper  design. 

The  maximum  line  tension  may  be  considered  as  either  the  combina¬ 
tion  of  dead  weight  and  dynamic  load  or  the  snap  load.  Since  the  dead 
weight  is  a  constant,  it  is  easily  determined.  The  maximum  dynamic 
load  is  defined  as  the  peak  of  the  tension  increase  or  decrease  from  the 
static  tension  load.  Dynamic  loads  are  generated  by  the  longitudinal 
and  horizontal  oscillations  at  the  surface  end  of  the  line.  The  snap 
load  is  an  impact  load  caused  b;  the  sudden  retensioning  of  the  line 
after  a  state  of  zero  tension  (Figure  1)  which  occurs  when  the  cable 
sy°tem  is  subjected  to  surface  motions  of  large  amplitude  and/or  high 
frequency.  Once  the  decrease  in  dynamic  tension  exceeds  the  magnitude 
of  the  static  load,  slack  occurs  in  the  line.  A  snap  load  develops  in 
the  line  as  the  platform  subsequently  moves  upward.  The  duration  of 
the  snap  load  is  short  but  its  amplitude  may  be  many  times  greater  than 
the  maximum  dynamic  load,  depending  upon  the  properties  of  the  line  and 
the  payload.  It  is  therefore  necessary  to  know  the  expected  peak  snap 
load  in  mooring  or  lifting  operations  in  order  to  provide  an  adequate 
cable  size. 

While  dynamic  tension  solutions  for  lifting  and  mooring  systems 
have  been  sought  by  many  investigators  (references  1-7)  the  snap  load 
problem  has  received  relatively  little  attention.  Goeller  studied  snap 
loads  for  a  single  degree  of  freedom  system  (reference  8).  Wilkins  de¬ 
veloped  hn  approximate  solution  for  the  problem  of  a  short  cable  system 
(reference  9).  Snap  tensions  caused  by  a  free  falling  anchor  have  been 
investigated  by  Nath;  good  results  were  achieved  using  a  lumped  muss 
model  und  solving  the  problem  with  a  predictor  -  corrector  type  numerical 
method  (reference  10).  However,  methods  for  predicting  snap  louds  in 
long  cable  systems  are  not  available,  nor  is  there  any  literature  on 
snap  loads  in  catenary  mooring  lines. 

The  presetu.  work  involves  the  development  of  a  computer  method  for 
the  calculation  of  snap  loads  in  lifting  und  mooring  lines  caused  by 
oscillatory  excitations  at  the  surface  (Figure  2).  Exact  solutions  of 


the  problem  are  not  feasible  because  a  cable  cannot  resist  compression. 
Therefore,  Hooke's  law  cannot  be  fulfilled  when  the  line  tension  drops 
below  zero.  This  nonlinear  elastic  property  of  the  cable  makes  the 
equation  of  motion  nonlinear.  In  this  work,  a  lumped  mass  model  is 
used  and  numerical  solutions  are  sought.  The  equations  of  motion  are 
written  for  each  mass  point  and  are  solved  by  a  computer  subroutine 
which  has  proven  to  be  effective  in  solving  a  set  of  nonlinear  first 
order  differential  equations. 

The  end  product  of  this  work  is  a  Fortran  IV  computer  program  called 
SNAPLG.  This  program  can  solve  transient,  dynamic  and  snapload  responses 
to  surface  excitations  of  any  composite  underwater  cable  system  in  any 
current  conditions  with  fixed  or  weighted  lower  end  boundary  conditions. 
The  program  is  applicable  to  deep  sea  emplacements  and  retrievals,  deep 
sea  mosjrings  and  the  construction  of  underwater  cable  systems.  It  is 
also  valid  for  analyzing  the  free  fall  anchor  dynamics.  The  limitations 
are  planar  motion  and  one  continuous  cable.  Three  dimensional  motion 
and  branched  out  legs  can  be  included  in  the  program  with  some 
modifications. 


NUMERICAL  SOLUTION 
Lumped  Mass  Model 

The  cable  system  is  divided  into  a  number  of  short  line  segments. 
Each  such  segment  is  assumed  to  be  small  enough  that  it  acts  as  a  rigid 
body,  whose  centroid  is  tailed  a  mass  point.  Furthermore,  all  external 
forces  act  through  this  point.  Adjacent  mass  points  are  connected  by 
springs  whose  elastic  property  is  equivalent  to  that  of  the  line  segment 
between  the  mass  points.  Dashpots  are  attached  to  mass  points  to  aim- 
ulate  hydrodynamic  damping  provided  by  the  fluid  drag  force  along  the 
line  segment. 

Any  cable  system  suspended  from  the  surface  may  be  modeled  as  shown 
in  Figure  3.  The  upper  end  represents  the  forced  circular  or  ellipsoidal 
surface  motions.  The  bottom  end  of  the  cable  may  be  a  heavy  concentrated 
payload  as  in  a  lift  system  or  an  anchor  fixed  to  the  seafloor  as  In  a 
mooring  system. 

The  mathematical  model  was  based  on  a  single  set  of  coordinates; 

■b  x  is  the  horizontal  distance  to  the  right  of  the  origin  and  +  z  is  the 
vertical  distance  above  the  origin.  The  basic  mass  point  geometrical 
relationashlp  is  shown  in  Figure  6.  The  origin  can  be  located  anywhere 
desired.  But  a  convenient  point  is  at  the  lower  end  of  the  unstretched 
cable. 
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The  distance  between  mass  points  and  the  inclinations  are  expressed 
in  terms  of  the  x  and  z  coordinates  as  follows: 


^i  “  V  ^Xi  xi-l^  +  ^zi  ~  Zi-1^2 
*i+i =  -  v2  +  (,i+i  -  v2 
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i 


arctan 
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arctan 
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arctan 


zi+l  ~  2i-l 
Xi+1  "  xi-l 


Mathematical  Formulation 

The  equation  of  motion  for  the  I*’*1  mass  is  derived  based  on  the 
principle  of  dynamic  force  equilibrium.  Therefore,  by  resolving  forces 
into  the  x  and  z  directions,  the  equilibrium  condition  for  the  ith  mass 
can  be  expressed  as  follows: 


<Vi +  <Vi +  <D,>1 +  <Vi«  +  <Vi  +  >s  - 0 


<Vi +  <Vi +  <D,>i +  <Vt«  *  <Vi +  F. 


(1) 


where  the  subscripts  x  and  z  denote  direction  of  the  force,  the  sub¬ 
scripts  i  and  i+1  identify  the  appropriate  mass  or  line  segment,  and 

I  -  inertial  force 

8  *  body  or  gravitational  force  (weight) 

D  *  hydrodynamic  drag  force 
T  *  cable  tension 
F  -  imposed  external  force 
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The  inertial  force  is  the  product  of  the  virtual  mass  and  the  accel¬ 
eration.  For  the  ith  mass  point  the  x  and  z  components  of  this  force  are 


(I  1  -  d2x 


.  d2I 

1  1  =  "i  ~2 
at 


(2) 


where 


“i  “  g  {°‘5  C1  li  +  0,5  '•i+1  *1+1  +  0-5  11  [jf,l  '  '■!  1  ‘l 


+  -  ‘W  h+l]  +  “l  +  *'  <Wi  -  wi  >} 

g  *  acceleration  due  to  gravity 
i  «  dry  weight  per  unit  length  of  cable 
C  -  submerged  weight  per  unit  length 
i°  «  length  of  segment  at  zero  tension 
a  ■  mass  coefficient  of  cable 

W  *  dry  weight  of  payload 
W'  “  submerged  weight  of  payload 

a'  *  mass  coefficient  of  payload 

c  »•  'itti! 


The  body  force  is  simply  the  gravitational  force  or  the  submerged 
weight  of  the  cable  segment  and  the  In-line  mass  at  the  i t5'  mass  point. 
Therefore,  the  x  and  a  components  are : 


^xh  “  0 


O) 


A 


The  current  drag  is  calculated  base^i  on  the  assumption  that  the 
vertical  component  of  the  current  velocity  is  negligible  and  that  the 
tangential  drag  along  the  cable  is  assumed  to  be  negligible.  The 
orientatio/  jf  t  le  payload  is  assumed  constant  with  change  of  cable 
slope.  The  relative  horizontal  velocity  of  the  mass  point  become^,: 


P 


i 


(A) 


There  are  three  hydrodynamic  damping  forces:  (1)  the  payload,  (2)  the 
normal  cable  drag  and  (3)  the  tangential  cable  drag. 


(Vt  •  -(cVi  Ml>i|-<Vi 810  n'Vtl'VihVi cus  VVibVj 


da . 

(t)  ) ,  -  -(C‘  >.  Ttr 
i  i  s  i  at 


(5) 


a.  r(v,  '«•  Wi  (Vi  -<ct>i sin  M<Vi  <Vi 


where  (C*  ).  -0.3  (s>  ‘  )  c  (A  ). 

A  1  ^  i  ■  *  4 


<cVi  *  #,J  (flVt  *'u  <Vt 


<Vi  - 0  “  ‘Vi  ■  «t  S  *  dt*t  W 


‘Vi  *  ‘Vi  •*„  *  wt  ‘i  *  Ji.l  ‘ft1 

d«. 


(V  ) 

'Vi 


“  pi siw  *i "  ic  m*  “t 
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<Vt  -  P4  co*  -  — -  sin 


-  drag  coefficient  of  cable 

£ ’  *  drag  coefficient  of  payload 

a  «  saws  density  of  seawater 
w 

A  -  payload  drag  area 
d  -  cable  diameter 
Q  *  current  velocity 


The  elastic  forces  are  calculated  based  on  Hooke's  Law  using  the  spring 
constant  of  the  cable  segment. 

(T  )  =■  -  k.  (2.  -  2.°)  cos  $ 

tor  (ti-  tl°)>o  (6) 

(Tz>i  a  -  (2^  -  4^°)  sin  ^ 

and 

Tt  =0  for  (2i  -  t*)  <  o 

Other  external  force  components  (F  and  (F„)^  are  to  be  defined  as 
known  quantities. 

By  substituting  aquations  (2)  through  (6)  into  the  corresponding 
portion  of  equation  (1),  the  equation  of  motion  of  the  ith  mass  point 
in  the  x  and  a  directions  can  be  obtained. 

Boundary  Conditions 

The  upper  end  boundary  eanditivns  may  be  expressed  as  follows? 


x  »  W .  ♦  r  sin  wc  +  o  e 
n  n  x  n 

a  *•  v  +  f  jjia 
ft  n  % 

wild  re  '  i,  e  *  steady  state  coordinates  in  s  and  z  direction* 

respectively 

t  ,  r  *  peaiv  upper  end.- displacement  in  s  and  a  direction* 
*  respectively 

v  *  excitation  frequency 

u^  *»  forced  upper  end  velocity  in  x  direction 

With  the  proper  cosbinatiott  of  values  of  end  uft,  surface  exci¬ 

tation  of  several  configurations  are  possible. 


The  lower  end  can  be  either  free  or  fixed.  Free  end  cable  systems 
may  have  a  heavy  payload,  a  buoy  or  no  payload  at  the  free  end.  For 
free  ends,  the  boundary  conditions  can  be  described  by  setting  i  =  o  in 
equations  (1)  through  (6).  The  resulting  equations  are: 


d2x 


dt‘ 


O  _1 

m 

o 


-  (C"  )  P  P  j  -  (CL)  sin  6  <V„)  O')  -  (C  ) 

x  o  o  oj  N  o  o  N  o  N  o  t  t 


.2 

d  z 


dt 


cos 

0 

(VJ 

(v  ) 

+ 

1 

0 

0 

t  ( 

dz 

3 

dz 

i 

(C'  ) 

Z  0 

o 

o 

m 

0 

dt 

dt 

kl  '*1 


X  o 


(CL)  cos  6  (VM)  (V  ) 
No  o  No  No 


(7) 


«Vo 


sin  0  (V  )  (V  )  +  k  (J sin  0  -0.55/4.  -  W'  +  (F  ) 

ototo  ill  c  11  o  zo 


where 


<Vo  -  °-25  ‘Vo  “l  h 
<Vo  -  °-25  <8t>.,  f„  » 

(CVo  *  °'5  «'x>o  <Vo 
(CVo  *  °'5  <6  Vo  ’*  ‘Vo 

<Vo  ■  po  8ln  Co  dT  008  eo 


dz 

(V  )  =>  P  cos  6  -y—  sin  0 

too  o  dt  o 

For  a  fixed  lower  end  at  the  origin  the  boundary  conditions  are: 


x  -  0 
o 


for  all  t's 


z  =0 
o 
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Initial  Conditions 


Physically  any  selected  initial  condition  will  eventually  lead  to 
a  steady  harmonic  response  to  the  external  excitation.  The  transient 
period  depends  on  the  frequency  response  of  the  system.  An  initial 
condition  which  describes  a  realistic  initial  state  of  the  cable  system 
requires  little  or  no  transient  period.  On  the  other  hand,  because  a 
numerical  solution  technique  is  used,  a  poorly  prepared  initial  condi¬ 
tion  can  result  in  instability  at  the  beginning  of  the  calcula-ion. 
Therefore,  it  is  most  desirable  to  determine  first  the  static  configur¬ 
ation  of  the  cable  system  in  question  under  static  loads.  Velocities 
of  the  mass  points  may  be  assigned  if  necessary  to  matci-  ,ie  upper 
boundary  condition.  Analytical  solutions  are  available  for  simple  lift 
and  mooring  systems  with  uniform  line  properties  but  with  no  in-line 
concentrated  masses  and  current  drag  forces.  Computer  programs  can  be 
used  to  calculate  the  static  state  configuration  of  complex  cable  sys¬ 
tem  in  a  current  field.  The  static  configuration  can  then  be  input  to 
Program  SNAPLG  as  the  initial  conditions  for  dynamic  analysis.  Although 
causing  a  somewhat  longer  transient  time  a  graphically  determined  static 
configuration  is  often  adequate. 


For  a  simple  lift  system,  the  static  configuration  is  a  vertically 
stretched  cable.  The  elongation  of  the  ith  line  segment  due  to  static 
loadings  of  system  dead  weights  and  other  static  forces  may  be  expressed 
as: 


‘i-1 

o 


k  [  +  ‘Vi  j 


The  coordinates  of  the  i£^  mass  point  become: 


zi  B  Zi-i  +  h  +  ei 


z  « 

o 


-*ei 


x 


i 


=  0 


For  a  catenary  shaped  cable  system  the  static  configuration  is  cal¬ 
culated  based  on  inextensible  catenary  equations.  Assuming  that  the 
upper  end  is  subjected  to  a  vertical  load,  V,  and  a  horizontal  load,  H, 
and  that  the  lower  end  is  fixed,  the  tension  at  the  upper  end  is: 

T2  =  H2  +  V2 
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Based  on  catenary  equations: 


H 


cosh 


T 

H 


and 


{'gn  -  V 

H 


T-H 

V 


where*  xh  and  zh,  are  the  coordinates  of  the  lowest  point  on  the 
catenary,  where  the  slope  is  zero;  xn  and  zh,  are  the  coordinates  of 
the  top  mass  point. 


H 


0 


and  the  lower  end  tension  is: 


T  =  -z,  £  '  +  H 
o  h 


therefore, 


the  total  scope  cf  the  cable  is  calculated  based  on  the  third  catenary 
equation: 


r 


sinh 


-  sinh 


The  imaginary  cable  scope  to  the  lowest  point 


of  the  catenary  is 


Let 


As 


s 

n 
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and 


=  sq  +  (i  -  0.5)As  for  i  =  1, 


.  n 


Sn+1  =  s 


Using  again  the  catenary  equations: 

H  -1  1  ‘l  - 

*i  *  r  3l“h  “  +  \ 


Z^  =  — T  (cosh 


-1  rcrl  -  V 


-  1)  +  z. 


for  i  ■  1, 


(n  +  1) 


Method  of  Solution 


There  are  two  differential  equations  for  each  mass  point.  For  n 
mass  points,  the  number  of  equations  is  2n.  These  equations  are  second 
order  second  degree  differential  equations  in  terms  of  x  and  z.  In 
order  to  solve  these  equations  by  the  Adams  numerical  method,  each  equa¬ 
tion  is  transformed  into  first  order  differential  equations  by  substituting 


v 


i 


2 

dt  dt 


dt2  dt 


As  a  result,  there  are  a  total  of  4n  first  order  parallel  differential 
equations  to  be  solved  by  numerical  integration. 

The  fourth  order  Adams-Predictor-Corrector  Method  which  is  a  method 
for  solving  a  system  of  first-order  ordinary  differential  equations 
(reference  12)  is  used  in  the  present  study.  The  Runge  Kutta  formula 
is  used  to  start  the  integration  and  whenever  the  time  step  size  is 
changed.  The  Adams-Bachforth  and  Adams-d-loulton  formulas  are  used  as 
predictor  and  corrector  respectively.  The  relative  and  absolute  errors 
are  computed  and  compared  with  the  predetermined  error  bounds.  If  the 
relative  and  absolute  errors  are  either  too  small  or  too  large,  the  time 
step  size  will  be  changed  automatically  to  reduce  computing  time  and 
instability.  It  also  has  a  provision  for  storing  and  printing  the  re¬ 
sults  at  selected  time  increments  (reference  13) . 
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A  snap  load  test  was  developed  to  detect  the  buckling  of  the  line 
segments  and  reassign  the  value  of  the  segment  spring  constant  to  ac¬ 
count  for  the  lack  of  compressive  stiffness.  The  test  is  made  before 
each  computation  of  the  velocity  and  acceleration  at  each  mass  point. 
The  line  buckling  is  detected  by  comparing  the  instantaneous  cable 
segment  lengths  with  the  unstretched  lengths.  If  slackness  is  detected 
in  a  segment,  the  corresponding  elastic  spring  constant  kf  is  set  to 
zero.  When  this  segment  becomes  taut  again  the  spring  constant  is  re¬ 
assigned  Its  original  value.  This  simulates  the  temporary  disassocia- 
tion  of  two  adjacent  mass  points  caused  by  a  slack  line  segment.  As 
these  two  points  move  far  enough  apart,  the  line  segment  regains  its 
ability  to  resist  tension.  The  basic  set  of  differential  equations  has 
thus  been  altered  to  fulfill  the  condition  of  line  buckling. 


PROGRAM  SNAPLG 

Computer  program  SNAPLG  consists  of  the  main  program  and  six  sub¬ 
routines.  These  subroutines  are  DIFEQ1,  ADAMS,  LOAD,  RNGKTA,  ERRTST  and 
STATIC. 

Program  SNAPLG  provides  entrance  for  line  and  payload  properties. 

It  also  calculates  the  initial  configuration  of  the  system  if  not  given 
as  input.  Properties  for  each  line  segment  to  be  input  are  segment 
length,  cross-sectional  area,  diameter,  dry  and  wet  unit  weights,  and 
Young's  modulus  of  elasticity.  The  properties  for  each  in-line  payload 
include  the  weight  in  air  and  in  water,  mass  and  drag  coefficients,  drag 
areas  in  vertical  and  horizontal  directions.  The  upper  boundary  condi¬ 
tions  are  the  vertical  and  horizontal  components  of  line  tension  at  the 
support.  Surface  displacement  excitation  is  input  in  the  form  of  ampli¬ 
tude  and  frequency.  Also  input  may  include  the  horizontal  velocity  of 
the  surface  ship.  The  depth  of  water  must  be  given  for  catenary  systems 
and  the  depth  of  suspension  must  be  given  for  lift  systems.  The  desired 
period  of  real  time  computation  must  be  specified.  The  current  velocity 
at  each  mass  point  may  be  included  if  required. 

Based  on  the  depth  of  suspension  the  unstretched  position  of  each 
mass  point  is  computed  for  a  lift  system.  The  positions  will  be  used 
as  the  initial  condition.  For  a  catenary  system,  the  initial  positions 
can  be  either  input  or  calculated  by  SUBROUTINE  STATIC  based  on  surface 
tension  and  line  properties. 

All  constants  used  in  the  differential  equations  are  computed  here. 
Complete  time  domain  solutions  can  be  printed  and  plotted  at  the  end  of 
the  SNAPLG. 
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SUBROUTINE  DIFGQ1  has  two  functions.  First,  it  tests  the  slack 
condition  in  each  line  segment  and  corrects  the  spring  constant  when 
necessary.  Secondly,  it  computes  the  top  boundary  conditions  and  cal¬ 
culates  the  acceleration  and  velocity  of  each  mass  point. 

SUBROUTINE  ADAMS'  main  function  is  to  select  the  best  time  incre¬ 
ment  for  the  next  round  of  calculations.  It  starts  the  calculations 
using  the  RUNGE-KUTTA  method  and  continues  the  calculations  with  the 
ADAMS  method.  It  also  stores  and  outputs  calculated  information  at 
desired  time  intervals. 

SUBROUTINE  LOAD  is  called  by  ADAMS  at  desired  time  intervals  when¬ 
ever  results  need  to  be  output.  The  output  includes  the  real  time, 
position  and  velocity  of  each  mass  point  and  the  line  tension  in  each 
segment.  This  information  is  also  stored  in  a  matrix  to  be  called  out 
when  needed  at  the  end  of  SNAPLG. 

SUBROUTINE  RNGKTA  is  called  by  ADAMS  to  calculate  the  position  of 
each  mass  point  at  the  beginning  of  the  calculation  and  after  each  out¬ 
put  time  increment. 

LOGICAL  FUNCTION  ERRTST  is  also  called  by  ADAMS.  It  determines 
whether  the  error  between  the  predicted  and  corrected  value  is  within 
the  desired  bound. 

SUBROUTINE  STATIC  is  called  by  PROGRAM  SNAPLG  to  compute  the  static 
configuration  of  a  catenary  system  based  on  the  assumption  that  the  line 
is  not  stretchable. 


RESULTS 

An  example  problem  was  analyzed  by  program  SNAPLG.  The  cable  was 
a  1-inch  wire  rope  mooring  line  suspended  in  600  feet  of  water.  The 
bottom  end  is  fixed  to  the  seafloor  and  the  top  end  is  subjected  to 
vertical  and  horizontal  sinusoidal  excitations.  Current  loads  were 
neglected.  Part  of  the  computer  output  is  plotted  in  Figure  7.  The 
snap  loads  are  the  sharp  tension  increase  immediately  following  each 
occurrence  of  slack.  The  second  peaks  are  the  dynamic  peaks.  In  this 
case,  the  value  of  peak  snap  loads  is  smaller  than  that  of  the  dynamic 
loads,  probably  due  to  the  large  damping  capability  of  the  catenary. 

The  snap  loads  are  expected  to  be  larger  for  larger  surface  excitations. 
The  output  is  reasonably  stable  and  consistent  at  the  end  of  the  second 
wave  period.  This  rapid  convergence  characteristic  means  a  savings  in 
computer  time. 

A  second  example  is  the  deployment  of  a  buoy/anchor  assembly.  The 
3800-pound  anchor  is  suspended  1800  feet  below  a  cylindrical  buoy  by  a 
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1-inch  Samson  Cordage  Power  Braid  rope.  The  buoy  has  a  net  buoyancy  of 
2500  pounds.  The  whole  assembly  is  supported  50  feet  under  the  platform 
by  a  5/ 8- inch  Samson  Cordage  Power  Braid  rope. 

The  platform  is  assumed  to  move  ±  2.0  feet  from  its  mean  position 
at  a  period  of  3.5  seconds.  The  result  is  plotted  in  Figure  8.  The 
three  curves  represent  the  total  tension  in  the  three  line  segmc  ,ts. 
Segment  12  is  the  5/8-inch  rope  above  the  buoy.  A  snap  load  occurs  only 
in  segment  12.  The  buoy  acts  as  a  filter  to  dampen  out  the  shocks  from 
the  surface.  The  values  of  the  snap  load  peaks  are  not  high  since  the 
synthetic  rope  can  absorb  a  great  amount  of  shock  energy. 

The  third  example  is  a  much  more  complex  cable  system  construction 
problem.  A  mooring  leg  is  being  deployed  in  deep  water.  The  leg 
consists  of  4186  feet  of  1/2-inch  diameter  3  x  19  three-conductor  torque 
balanced  wire  rope  which  weighs  0.49  Ib/ft  dry  and  0.31  lb/ft  wet.  The 
cable  is  supported  by  a  buoy  weighing  525  lb  dry  and  has  a  net  buoyancy 
of  1625  founds.  Attached  at  the  lower  end  is  a  25,000  pound  clump  anchor 
weighing  20,000  pounds  submerged.  A  crown  line,  25,000  feet  of  1-inch 
3  x  19  three-conductor  wire  rope,  is  used  to  lift  and  position  the  anchor 
to  the  desired  site.  It  weighs  1.75  Ib/ft  in  air  and  1.23  lb/ft  in 
water.  The  upper  end  of  the  crown  line  is  supported  by  a  surface  ship 
which  is  subjected  to  wave  actions.  In  this  example,  the  ship  motion 
consist?  of  a  constant  velocity  of  3  knots  pulling  away  from  the  buoy 
and  a  circular  motion  of  3  foot  radius  and  1.26  radians  per  second  fre¬ 
quency  in  a  vertical  plane.  A  horizontal  force  of  810  pounds  is  exerted 
on  the  buoy  by  the  balance  of  the  cable  system.  The  dynamic  effect 
caused  by  the  surfaced  motion  is  assumed  negligible  at  the  buoy.  Figure 
9  is  a  schematic  presentation  of  the  cable  system  in  the  problem.  Figures 
IQ-17  show  the  tension  history  of  each  segment  in  the  system.  Although 
the  displacements  and  velocities  in  each  mass  point  are  not  shown  here, 
they  can  be  displayed  in  the  same  fashion  as  the  tensions. 

DISCUSSION 

Program  SNAPLG  is  a  useful  tool  for  designing  hardware  and  planning 
deep  ocean  deployment  and  recovery.  The  snap  load  is  generally  affected 
by  four  major  elements:  Seastate,  surface  platform,  cable  and  payload, 
or  static  tension.  With  these  elements  as  input  to  the  program,  the 
occurrence  of  snap  load  may  be  predicted  and  the  magnitude  of  the  snap 
load  caluclated.  If  based  on  the  maximum  tension,  the  cable  is  not  safe, 
any  of  the  four  inputs  may  be  modified  for  another  trial.  The  iteration 
can  continue  until  the  most  favorable  combination  of  four  elemencs  has 
been  reached.  For  a  lift  system,  the  maximum  tension  should  be  checked 
throughout  the  depth  of  the  emplacement  and  recovery.  Particularly, 
the  depth  where  a  prolonged  payload  suspension  is  anticipated.  For  a 
catenary  mooring  cable,  the  tension  for  a  series  of  catenary  configura¬ 
tions  should  be  checked. 
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The  surface  excitation  is  input  as  sinusoidal  motions  of  a  surface 
platform.  By  modifying  the  equations  of  motion,  random  motion  may  be 
input  as  a  Fourier  series.  The  near  surface  wave  particle  velocity 
loadings  may  also  be  included  into  the  equation  of  motion  without 
difficulty.  This  method  may  be  useful  in  solving  the  cable  strumming 
problem.  The  only  difficulties  will  be  the  determination  of  the  current 
induced  lift  force  on  a  strumming  cable  and  tee  large  number  of  mass 
points  required. 

Information  needed  on  the  line  properties  and  the  payload  proper¬ 
ties  are  standard.  It  is  desirable  to  know  the  static  configuration  of 
a  catenary.  If  this  is  not  available,  the  top  tension  and  the  inclina¬ 
tion  should  be  given.  An  approximate  static  configuration  can  be  ob¬ 
tained  by  a  graphic  method. 

By  using  proper  control  parameters,  the  program  will  select  the 
proper  subroutines  for  the  lift  or  moor  calculations.  Control  variables 
for  the  ADAMS  subroutines  include  tee  initial  time  step  size,  relative 
and  absolute  error  bounds,  length  of  real  time  calculation,  and  the  time 
intervax  of  printing  output.  The  results  of  the  computation  are  plotted 
on  oscillographs  for  presentation  and  visualization.  Movies  can  be  made 
from  these  graphs  to  animate  the  dynamics  of  the  cable  system. 

One  advantage  of  this  program  is  that  the  compute”  versus  real  time 
ratio  is  less  than  2.  This  is  based  on  a  limited  number  of  runs  on  a 
5  segment  cable  system.  A  twenty-second  real  time  computation  can  be 
made  with  a  cost  of  no  more  than  $40.00  using  a  CDC  6600  computer. 

This  method  can  be  extended  to  three  dimensional  cable  systems  with 
branched  out  legs  and  members.  Three  differential  equations  will  be 
required  Jo  define  the  motion  of  the  mass  point.  Consequently,  longer 
computation  time  is  expected.  There  is  no  forseen  technical  difficulty 
in  achieving  the  three  dimensional  capability. 

The  accuracy  of  the  solution  depends  on  the  mathematical  model 
adopted  for  the  problem.  High  precision  may  be  achieved  with  a  large 
number  of  lumped  masses,  small  error  bounds  and  small  output  time 
intervals.  This  program  has  not  been  validated  by  experimental  data. 

CONCLUSIONS 

1.  A  computer  program  SNAPLG  ha1)  been  developed  to  solve  the  maximum 
peak  tension  in  a  lift  line  or  in  a  mooring  line  subjected  to  surface 
excitations.  The  program  is  based  upon  a  numerical  method  called  the 
Fourth  Order  ADAMS  Predictor-Corrector  Method.  Lumped  mass  models  were 
used  to  simulate  the  cables.  Hie  program  tests  the  line  tension  contin¬ 
uously  to  avoid  negative  values.  Consequently,  snap  load  conditions  can 
be  simulated  and  predicted. 
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2.  The  program  provides  a  useful  tool  for  the  design  of  hardware  and 
the  emplacement  and  recovery  of  heavy  payloads  and  deep  mooring  systems. 
Knowledge  of  the  expected  dynamic  and  snap  loads  makes  possible  the 
optimal  design  of  the  lift  and  mooring  systems  and  the  handling  equipment. 
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x  =  horizontal  coordinate  (ft) 

x  =  horizontal  coordinate  under  steady  state  conditions  (ft) 
z  =  vertical  coordinate  (ft) 

2  *  vertical  coordinate  under  steady  state  conditions  (ft) 
a  3  mass  coefficient  of  cable 

=>  mass  coefficient  of  payload 

3  »  drag  coefficient  of  cable 
8'  *  drag  coefficient  of  payload 
i  ■  unit  weight  of  cable  (lb/ft) 

£■*  »  submerged  unit  weight  of  cable  (Ib/ft) 

t  *  elongation  of  cable  (ft) 
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Subscripts: 

h  *  refers  to  location  on  catenary  where  slope  Is  zero 
l  *  refers  to  I6*'  mass  point  or  i£*‘  cable 
N  “  component  normal  to  the  cable 
n  *  upper  end  of  cable 
o  *  lower  end  of  cable 
t  “  compoaeut  tangential  to  the  cable 
x  •  component  in  x  direction 


«*  component  in  z  direction 


Figure  i,  Static  toad  i;>  the  weight  of  the  cable  system 
Dynamic  load  i.i  additional  load  caused  by  the 
.motion  of  the  system.  Snap  load  is  usually 
the  maximum  possible  load  in  the  system.  It 
is  created  by  a  sudden  tensioning  of  a  slack 
cable  system. 


Figure  ?,  Configurations  of  lift  and  see r  cable  systeus 
subject  to  surface  excitations. 


Figure  5 


z  i+1  (elastic  force) 


Figure  6.  Static  configuration  of  a  catenary  cable  system. 
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Figure  7.  Snap  load  in  a  1-inch  diameter  wire  rope  mooring  line  1400  ft.  long. 
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Figure  8.  s„,p  load  in  a  buoy-„chor 
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